Stream food webs in tropical mountains rely on allochthonous carbon regardless of land use

The relative importance of allochthonous and autochthonous carbon (C) as sources of energy for tropical stream food webs remains an open question. Allochthonous C might be the main energy source for small and shaded forest streams, while autochthonous C is more likely to fuel food webs draining land uses with less dense vegetation. We studied food webs in cloud forest streams draining watersheds with forests, coffee plantations, and pastures. Our goal was to assess the effects of those land uses on the C source and structure of stream food webs. The study took place in tropical montane streams in La Antigua Watershed, in eastern Mexico. We selected three streams per land use and sampled biofilm and leaf litter as the main food resources, and macroinvertebrates and aquatic vertebrates from different trophic guilds. Samples were analyzed for δ13C and δ15N isotopes. Using a Bayesian mixing model, we estimated the proportional assimilation of autochthonous and allochthonous carbon by each guild. We found that consumers were mostly using allochthonous C in all streams, regardless of watershed land use. Our findings indicate that montane cloud forest streams are dominated by allochthony even in watersheds dominated by pastures. Abundant precipitation in this life zone might facilitate the movement of allochthonous C into streams. While food webs of streams from coffee plantations and pastures also rely on allochthonous resources, other impacts do result in important changes in stream functioning.


Introduction
Stream food webs rely on two major sources of C, allochthonous inputs from terrestrial vegetation (e.g., leaf litter) and autochthonous inputs of algal growth in the channel [1].The type and amount of riparian vegetation and catchment land use determines the amount and quality of resources available for stream food webs.In most temperate regions, inputs of riparian C are the main energy source for small forest stream food webs [2].Shading from riparian vegetation limits algal growth, reducing inputs of autochthonous C and multiple trophic levels are linked to allochthonous C [3].Studies in tropical streams are less definitive, with some authors arguing for a larger role of algal C as the main source of energy for tropical food webs [e.g., 4,5].Studies combining measures of food consumption with stream metabolism indicate that the overall pattern could be similar to that described for temperate regions [6].Interestingly, studies in small tropical streams in Brazil indicate that the amount of primary productivity in forested streams is high enough to fuel aquatic food webs, but most consumers rely on allochthonous C [7].
Removal of riparian or catchment vegetation changes the proportions of allochthonous and autochthonous C available to consumers, severely altering stream food webs.Land uses that decrease canopy cover and reduce shading over the stream (e.g., pastures, logging) could increase algal production and decrease vertical inputs of leaf litter [8].Stream food webs respond to those changes by switching consumption from allochthonous to autochthonous C.This response has been documented in temperate regions [9] and a similar trend has been suggested for tropical regions, at least for some consumers (benthic macroinvertebrates [10]).Extreme land use changes (e.g., conversion of forest into pasture) cause the largest changes in stream food webs.Other land use changes (e.g., those that retain wide strips of riparian forest) have less dramatic effects.In tropical regions, less intense agricultural activities often replace riparian vegetation with other shrubs.Coffee plantations, for example, tend to leave no native vegetation at the stream edge, but coffee plants (and shade trees when present) shade stream channels and provide leaf litter for consumers [11].
Tropical montane cloud forest is a unique and valuable ecosystem that is changing rapidly due to human activities and land use conversion [12].Cloud forests are restricted to the upper parts of tropical mountain ranges, supporting high levels of endemism [13,14] and rich hydrological resources [15].They are a conservation priority for the protection of biodiversity and their rich hydrological resources [14,16].Land use conversion, from forest to agriculture (e.g., coffee) and pastures, along with a changing climate that might reduce cloud cover, are major anthropogenic impacts affecting this region [15,17].Cloud forests provide a unique opportunity to understand how land use changes affect stream food webs by altering their energy basis.Forest streams are often heavily shaded by riparian vegetation and experience a low light environment, with frequent cloud cover further limiting light penetration [18].Conversion into pastures drastically changes the light environment, potentially increasing algal C availability, while other land uses could create less drastic changes.Shade coffee plantations often reach stream margins and maintain low light levels, but impact streams by increasing erosion and consequently sediment, nutrient, and agrochemical inputs [18].
In this study, we focused on cloud forest streams draining watersheds with three dominant land uses that are common in many montane tropical regions: forests, coffee plantations, and pastures.Our goal was to assess the effects of those land uses on aquatic food webs and their main C sources in tropical montane streams in the La Antigua Watershed, Mexico.We hypothesized that the main C source is driven by the availability of autochthonous C. We predicted that forest streams would rely on allochthonous inputs of C as predicted for small forest headwater streams in general.Pasture streams were expected to have high algal growth and food webs based on autochthonous C. Coffee plantation streams should fall in between, but would be expected to be more similar to forest streams, because they are shaded and structurally alike to forest.

Ethics
Our study did not involve endangered or protected species.Vertebrate sampling was conducted under sampling permit FAUT-0303 from Secretarı ´a de Medio Ambiente y Recursos Naturales (SEMARNAT), Mexico.

Study sites
Our study was conducted in streams draining Tropical Mountain Cloud Forest (TMCF) (between 1500 and 2000 m a.s.l.) in the La Antigua watershed (Fig 1).La Antigua is a 2,623 km 2 basin that drains part of the states of Puebla and Veracruz (19˚050' and 19˚340'N; 96˚060' and 97˚160'W) in eastern Mexico, discharging into the Gulf of Mexico [19].Temperature and precipitation vary with elevation; annual means are 22˚C and 1300 mm at 600 m a.s.l. and 12˚C and 2000 mm at 2500 m a.s.l.[20].This watershed is an area of high biodiversity and hydrological importance [21].TMCF is the most biodiverse life zone in Mexico [22,23].Land use comprises a mixture of forest (mostly restricted to steep slopes), agricultural fields, shade coffee plantations, pastures, and urban areas [17].
We selected nine first or second order streams that are part of our ongoing research program in TMCF.At each stream, we delimited a 100-m reach that was dominated by one of the three study land uses: forests, coffee plantations, and pastures (Fig 1).Watersheds contained a mixture of land uses, but our selected reaches had at least a 100-m wide strip of the particular land use on each side.Forest streams had dense riparian vegetation of either natural forest or secondary growth.Coffee plantations were interspersed with trees (e.g., shaded coffee plantations), with the plantation often reaching the stream margin.Pastures were for cattle grazing and some riparian vegetation was present (e.g., scattered riparian trees).Stream physical and chemical characterization and aquatic organisms sampling were conducted within the 100 m study reach (Table 1).

Environmental variables and chlorophyll
Stream physicochemistry was measured at each end of the study reach following methods in Va ´zquez et al. [24] during the dry season of 2019.Mean canopy cover over the stream was estimated with the HabitApp mobile application by averaging canopy coverage of 10 photographs taken with a leveled cell phone at points separated ~10 m along each study reach.The channel and riparian zone slopes were measured using a clinometer.Temperature, dissolved oxygen, and electrical conductivity were determined in situ using an YSI multi-meter (Mod.85), and a potentiometer (Barnant Mod.20) was used to measure pH.Water samples were collected in 1 L glass bottles at each end of the study reach.Samples were transported to the laboratory at 4˚C and filtered through Whatman GF/C filters for analysis.In the laboratory, the following variables were measured using spectrophotometric techniques [25]: total suspended solids (TSS, dry weight), alkalinity (Alk, as CaCO 3 , phenolphthalein), ammonium (N-NH 4 + , Nessler), nitrate (N-NO 3 -, brucine), nitrite (N-NO 2 -, diazotization), chlorides (Cl-, titration with AgNO 3 ), sulfates (SO 4 , turbidimetric technique) and silica (SiO 2 , molybdate).Samples were collected in 250 ml glass bottles for measuring total phosphorus (TP, persulfate digestion, and ascorbic acid) and reactive phosphorus (RP, ascorbic acid).Calcium (Ca 2+ ) and magnesium (Mg 2+ ) were measured using an atomic absorption spectrophotometer (Varian Mod.AA240FS) and sodium (Na + ) and potassium (K + ) were measured with a flame photometer (Corning Mod.410).
Chlorophyll a concentration was measured at each study reach by collecting three rocks covered with biofilm and placing them in 90% methanol.Rocks were kept in the dark and refrigerated for 24 h before analysis.The chlorophyll a concentration was measured spectrophotometrically and its concentration was expressed in mg per m -2 [26].

Resources and consumers
Resources consisted of terrestrial plant tissue and aquatic biofilms.Plant tissue consisted of freshly fallen senescent leaves from outside the stream, riparian trees (i.e., green leaves), fern leaves, and pasture leaves.Samples were labeled and transported to the laboratory in plastic bags for processing.Biofilm was scraped from the exposed side of four or five rocks, using a small brush and tray.We collected two to three samples of biofilm per stream.Scraped material was moved into glass containers and transported to the laboratory.Consumers consisted of 10 aquatic taxa from all major trophic guilds that were found in most streams (Table 2).Organisms were collected with aid of kick nets by sampling the stream substrate along the study reach.All organisms were placed in containers with water, maintained in cool conditions, and transported alive to the laboratory.Trophic guilds were assigned following Merritt et al. [27].
All samples were transported to the laboratory for processing and drying at 70˚C for 48 hours.Leaves were placed in the oven in labeled paper bags, biofilm samples were left in their glass jars and oven-dried.Consumers were kept alive for 24 h at room temperature to allow for gut clearance, and were sorted by group and oven-dried for 48 h.Samples were ground to a fine powder, either manually or using a ball mill (Spex 8000D Mixer/Mill), for microcapsule preparation.All C and N isotope analyses were conducted at the Center for Applied Isotope Studies, the University of Georgia.Ratios of 13 C/ 12 C and 15 N/ 14 N are expressed in δ notation as the relative difference (‰) between the samples and international standards (Pee Dee Belemnite for C and air for N).
Resources and consumers were sampled during the dry and wet seasons of 2019.We did not assess seasonality, but expected a limited effect on isotopic signals.Seasonality in our study region affects the amount of resources available in streams, but not the types of resources [18] or composition of macroinvertebrate assemblages [28].

Stable isotope mixing model
The proportion of terrestrial organic matter (allochthony) on consumers was estimated individually for each taxa belonging to four consumer guilds (i.e., herbivore, detritivore, omnivore, and carnivore) using a system of two linked Bayesian models.The first model was used for the determination of algal stable isotope ratios.The second model was used to determine isotopic mixing.
In our model to determine algal isotope ratios, we estimated the δ13C and δ15N of the algal component of biofilm by removing the contribution of terrestrial sources (e.g., leaf litter).While various techniques have been proposed to isolate algae in biofilms [29], our approach has the goal of determining the degree of allochthony in the system, not the specific algal isotope signals.To delimit the possible range in algal values, we used as informative priors the mean and standard deviation of 35 previous measurements of algae in tropical streams (S1 Table in S1 Appendix).Published algal information is only used to delimit possible algal isotopic values.Similar approaches have been found informative in studies of lake food webs [30][31][32].
We estimated isotopic ratios of algae (δ 13 C and δ 15 N) and corrected isotopic ratios of leaf litter fitting a mixing model for two sources within a hierarchical Bayesian framework [33,34]: where δ P is biofilm stable isotopic signature (i.e.δ 13 C or δ 15 N), δ A is the algal stable isotopic signature, and F T is a terrestrial fraction on biofilm, and δ T is for stable isotope signature of allochthonous matter (i.e., leaf litter).We were most interested in the values of δ A for each land use; given that there were 3 categories (f, c, and p, for forest, coffee plantation, and pasture, respectively) [34].We set forest category as baseline (δ f ), such that its prior probability was where δ f for δ 13 C had a prior distribution of δ 13 C~Ɲ(−23.72,4.10), and for δ 15 N had a prior distribution of δ 15 N~Ɲ(4.10,3.13).Terms β c and β p correspond to group indicators for coffee plantations (c) and pasture (p), respectively.Terms for each land use category other than forest had a prior normalized distribution with a mean of 0 and a standard deviation (sd) of 1.As with δ A , we set forest category as a baseline for F T and used α c and α p as terms for coffee plantations and pasture respectively, such that where F f had a prior distribution of F f ~U(0,1), where 0 would be 100% of the autochthonous composition and 1 a 100% of the allochthonous composition.Terms α c and α p had a prior normalized distribution with a mean of 0 and a sd of 1.
The second model was a mixing model.After estimating the isotopic ratio of algae using the model described above, these values were passed into a mixing model of the form: where δ t,c is the isotopic ratio of each freshwater consumer (c) for each tracer (t), δ t,1:s is the vector of isotopic ratios for each of its potential carbon sources (1:s), Δδ is the trophic fractionation factor by each tracer (i.e., Δδ 13 C and Δδ 15 N), L is the trophic level for each consumer, and F is a vector of source proportions (drawn from a Dirichlet distribution with alpha set to a vector of ones).The symbol � represents the scalar product.δ t, 1:s had a prior normalized distribution with mean and sd values of δ 13 C or δ 15 N of our leaf litter samples, and estimated mean and sd for algae by each land use (S2  [35].L had a uniform distribution with prior probability between 0 and 10.We fit these models by simulating the posterior parameter distributions using the program Stan [36] using the rstan package in R [37].Stan uses a Markov chain Monte Carlo (MCMC) method to simulate posterior distributions.MCMC sampling of posteriors was performed with four chains and 5000 iterations by chain, after burn-in.We visually checked the chains for convergence and that of the scale reduction factor, R < 1:1, for all parameters.

Statistical analyses
Tendencies of variation in physicochemical variables among streams were assessed using a Principal Component Analysis (PCA).Variables were log 10 (x + 1) transformed to normalize the data, except for pH [38].The PCA used centered and standardized data.Analyzes were carried out with the MVSP program [39].Depending on the normality and homoscedasticity of data, we performed ANOVAS or Kruskal-Wallis One Way Analysis of Variance on Ranks to assess physicochemical differences among streams.

Land use and stream physicochemical characteristics
Stream physicochemistry during the dry season rendered few differences between streams (Table 1).There were no significant differences among land uses in most of their physical and chemical characteristics.Mean water temperature in coffee plantations was higher (17.1˚C) than in forest and pasture streams (14.6˚C).All streams were well-oxygenated (> 5.8 mg/L).Water conductivity in forest was higher than that of coffee and pasture streams (Table 1).The pH fluctuated from 6.6 to 7.0 in all streams.Significantly higher concentrations of Total Suspended Solids (TSS) were recorded in coffee plantation streams than in forests (H = 6.5, P <0.03).N-NH 4 concentrations were significantly higher in forests and coffee plantations (2.7-3.3 μM) than in pastures (0.1 μM) (H = 6.3,P <0.03).The concentration of N-NO 3 was higher in the forests than in the coffee plantations and pastures.The TP was recorded in higher concentrations in coffee plantations than in forests and pastures (Table 1).Chlorophyll a concentration was variable, but not significantly different among land uses (Table 1).
Streams draining each land use had similar physicochemical characteristics (Fig 2).PCA analysis grouped all streams by land use, except a single stream (Fig 2).Axes 1 and 2 explained 42% and 26% of the total variance.Forest streams and one of our pasture streams grouped in the quadrants on the right side of the biplot, displayed the highest conductivities and concentrations of N-NO 3 , oxygen, pH, SiO 2 , and chlorophyll a.In contrast, the other two pasture streams grouped in the quadrants on the left side of the biplot, showed low concentrations of all variables.The coffee plantation streams were grouped in the upper quadrants mainly associated with high TP, TSS, and temperature concentrations.

Basal resources and consumers
Biofilm and leaf litter were the basal resources sampled at all streams.Biofilm had similar isotopic values for C in forest than in pasture, but was highest in coffee plantations, while N was highest in coffee plantations.Leaf litter had similar C signals in forest and coffee plantations, and lower in pastures, while N was lowest in coffee plantations (Table 2).The isotopic position of biofilm was always in between leaf litter and algae, as it combines both resources (Fig 3).To address our main objective of determining the degree of allochthony on stream food webs, we focused on algal vs. leaf litter C isotopic signatures.We found little overlap between algae and leaf litter C signals in forest and coffee plantation streams; with pasture streams showing a larger degree of C overlap (S3 Fig in S1 Appendix).
Not all consumer taxa were abundant in all streams, but we collected enough samples of all trophic guilds (Table 2).Herbivores had similar C and N signals among land uses, average C ranged from -26.9 to -27.7‰ and N ranged from 4.3 to 4.6‰.Detritivores also had similar C among land uses (range: -26.0 to -28.3‰), but their N was lowest in pastures (range: 0 to 3.3‰).Omnivores also had a narrow range in their C signals (range: -25.5 to -27.9‰) and their N was variable with a larger range in pasture streams (range in pasture: 1.8 to 7.4‰).Carnivores had C signals that ranged from -24.6 to -26.5‰, and their lowest N signals occurred in pasture (range: 2.0 to 2.9‰; Table 2).
Bi-plots for stable isotopes, after accounting for trophic fractionation, showed that all trophic guilds had isotopic signals similar to leaf litter rather than algae in all land uses (Fig 3).As could be expected, detritivores are located close to leaf litter, while omnivores and carnivores are closer to biofilm in all land uses.Overall, consumers fall within the region covered by resource variability (Fig 3).
All trophic guilds show high levels of allochthony (Table 3).Consumers in streams draining the three land uses assimilated on average 34% of their C from leaf litter relative to only 20% from algae (means in Table 3).The degree of allochthony on stream food webs was determined by comparing assimilation patterns of algal C vs. leaf litter C. Detritivorous consumers showed the expected pattern of relying mostly on leaf litter and only decapods show some degree of algal consumption in forest streams (Fig 4).Herbivores in forest had low levels of leaf litter assimilation, but leaf litter was their main resource in coffee plantations and pastures (Fig 5).Omnivores consumed mostly leaf litter or had similar proportions of algae and leaf litter (Fig 6).Carnivores showed a reliance on leaf litter C, with some taxa also using algal C as an important resource (e.g., Megaloptera; Fig 7).

Discussion
Cloud forest streams in La Antigua watershed have food webs that rely primarily on allochthonous C, regardless of the land use surrounding them.Our findings indicate that food webs in forest, coffee plantations, and pasture streams have allochthonous and autochthonous resources available, but consumers mostly use C derived from allochthonous sources.This contradicts our initial hypothesis that stated that only forest streams would rely on allochthonous C, while pasture streams were expected to use mostly autochthonous C. Remaining patches of vegetation in riparian zones and high runoff in these rainy ecosystems potentially explain the persistence of terrestrial C. Our findings indicate that cloud forest streams have high levels of allochthony, based on the assimilation of C from terrestrial sources (following [30]).
The amount of allochthony in stream food webs is expected to be the result of land use cover and stream size [40].Small headwater streams draining forest receive ample inputs of leaf litter and their narrow channels do not create gaps in the canopy, limiting the amount of sunlight that reaches the stream bottom [40,41].The role of allochthonous C in fueling food webs in low order streams is well-established, in particular for temperate region streams [2,3].Multiple trophic levels in those streams are dependent on inputs of leaf litter from riparian vegetation [2].In contrast, studies in tropical streams surrounded by forest have produced mixed evidence [5].While some have documented similar degrees of allochthony as in temperate regions [42], others propose a larger role for algal C [4,43].Recent studies using a combination of approaches (e.g., stable isotopes and stream metabolism) has shown that some tropical stream food webs rely mostly on allochthonous C, even when algal C is available [7].
Our study provides further evidence for the importance of allochthonous C in cloud forest streams, where food webs are strongly linked to C from riparian zones.Streams draining natural areas of tropical cloud forest can be expected to rely on allochthonous resources.Our previous studies in this cloud forest area provide evidence for a high degree of allochthony in stream food webs.Using gut content analysis of dominant aquatic insects, Garcı ´a et al. [16] reported a predominance of plant tissue and amorphous detritus in most benthic consumers.The study found a range of 31-99% of aquatic insect gut contents composed of allochthonous material in forest streams.Similarly, gut contents in a pasture stream without riparian vegetation were found to have up to 82% allochthonous material.Gut content analysis has the limitation of focusing on ingested, rather than assimilated energy.However, our analysis using isotopes support those findings, suggesting that allochthony is dominant in forest and pasture streams in cloud forest.Further evidence for the importance of allochthonous resources comes from insect scrappers.Scrappers inhabit hard substrates in streams and are expected to consume mostly diatoms.However, this functional group was found to be a dominant component of macroinvertebrate assemblages in cloud forest streams https://doi.org/10.1371/journal.pone.0295738.g004[28].Their gut contents often contained high proportions of allochthonous material in both forest and pasture streams [16].
Deforestation of riparian vegetation and land use conversion from forest to agriculture alter many environmental characteristics and the functioning of stream ecosystems [8].Removal of riparian vegetation often changes resource availability for stream food webs, decreasing inputs of leaf litter by replacing riparian trees with crops or pastures [44].At the same time, a lack of arboreal vegetation increases solar radiation in the stream channel and enhances algal biomass and production [45].Macroinvertebrate assemblages often change their composition; losing groups specialized in consuming leaf litter and gaining scrappers and herbivorous taxa.This pattern has been described for temperate [46,47] and tropical streams [48,49].Our study shows that consumers continue relying on allochthonous C, even when coffee plantations and pastures provide reduced riparian cover relative to a forest.In our cloud forest study area, forest conversion to agriculture and pasture is a dominant trend [19] with consequent negative impacts on water physicochemistry [18] and macroinvertebrate assemblage composition and structure [28].

Land use and cloud forest streams
Land use had a clear effect on the physical and chemical characteristics of our study cloud forest streams.Our sampling of stream physicochemistry reflected conditions under reduced flow and few floods.Important differences were evident in the total suspended solids and ammonia concentrations in streams, which were more abundant in coffee plantation or pasture streams than in forest streams.Similar results have been reported in field studies and simulated analysis using SWAT models showing that land use increases TSS and nutrient concentrations in streams surrounded by agricultural lands compared to forests [50,51].Our results evidenced weaker land use impacts on algal biomass than those reported for the same area in previous studies [18].Algal biomass (i.e., chlorophyll concentrations) was similar among streams in our study, even slightly higher in forest.However, Va ´zquez et al. [18] found consistently higher algal biomass in pasture streams than in forest.The difference could be related to the much higher concentration of suspended solids that we found in our study.
The predominance of food web allochthony under all land uses suggests that leaf litter might not be the only source of terrestrial C for streams.Forest streams receive abundant direct inputs of C from riparian vegetation, via vertical inputs of leaves [28].Coffee plantations maintain riparian cover, even when it might be reduced or replaced by crop plants, and leaves inputs could be expected to be abundant.Our pasture streams maintained patches of riparian vegetation, which included scattered and isolated riparian trees and thin strips of riparian forest.However, in this wet and rainy environment, water runoff into the stream channel (superficial and as local groundwater) is abundant, in particular during the wet season [18], and streams can receive large amounts of fine particles from the stream banks and adjacent watershed.Allochthonous C could enter streams as fine particles during rain events, thus supplying food webs with organic matter.Scattered tree cover is often present in pastures, either close to the stream or as shade for cattle.This vegetation along with pastures could be an important source of C that is captured by the biofilm in streams and moves into food webs as scrappers and other consumers feed on the biofilms over rocks.
Our study documented, for the first time, the impact of land use on the C basis of stream food webs in a tropical cloud forest.The predominance of allochthony suggests that deforested watersheds maintain enough sources of allochthonous C to fuel stream food webs.The abundant precipitation is expected to be responsible for moving allochthonous C to the stream, even in open pasture areas, where stream consumers use it.In our study region, conversion of forest into other land uses results in significant changes in stream water physicochemistry, in particular increasing sediment and nutrient inputs into streams (our study; [18]).Biotic assemblages, especially macroinvertebrates, become simplified, as they lose sensitive taxa [28,52].Ecosystem processes are also affected by the land use surrounding the stream.For example, leaf-litter decomposition rates increase as forest cover decreases, potentially the result of increases in nutrient levels and other disturbances [53].Thus, while food webs remain based on allochthonous resources, multiple other ecosystem components change with land use change.

Fig 1 .
Fig 1. Sampling sites in the upper La Antigua watershed, Veracruz, Mexico.The nine study streams were first or second order draining forest (green stars), coffee plantations (red pentagons), and pastures (yellow circles).Since we selected first to second order streams, not all symbols fall on a visible stream channel.Figure modified from Va ´zquez et al. (2023).https://doi.org/10.1371/journal.pone.0295738.g001

Fig 2 .
Fig 2. Principal Component Analysis (PCA) of study sites based on water physicochemistry.Variables were measured at each end of the study reach and two points are representing each stream (points a and b).Numbers correspond to stream names in Table2and symbols to land uses, forest (green stars), coffee plantation (red pentagons), and pastures (yellow circles).https://doi.org/10.1371/journal.pone.0295738.g002

Fig 4 .
Fig 4. Resource consumption by detritivores.Posterior probability density of proportional source contributions by land use.The unshaded region corresponds to the 95% credible interval.

Table 1 . Physicochemical characteristics of the study streams, grouped by land use. Values
are means, standard deviations (SD), and minimum-maximum (Min-Max) values for the three streams sampled per land use type.Lines in bold represent significant differences among land uses, means with different letters are significantly different from the others.